THE MASS PROFILE OF THE COMA GALAXY CLUSTER 
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ABSTRACT 

We use a new redshift survey complete to m# ~ 15.4 within 4.25° from the center 



of the Coma cluster to measure the mass profile of the cluster to r ~ 5.5/i 1 Mpc. We 
extend the profile to r ~ 10ft, Mpc with a further sample complete to m/j = 15.4 in 
42% of the area within a 10° radius and to mz w = 15-5 in the remaining area. Galaxies 
within this region are falling onto the cluster on moderately radial orbits and thus 
do not satisfy virial equilibrium. Nonetheless, identification of the caustics in redshift 
space provides an estimate of the gravitational potential at radius r and hence of the 
system mass, M(< 10/i _1 Mpc) = (1.65±0.41) x 10 15 h~ 1 M & (1-a error). Previous mass 
estimates derived from optical and X-ray observations are limited to r < 2.5ft -1 Mpc. 
(T) ■ Our mass profile is consistent with these estimates but extends to distances four times 

as large. Over the entire range, the mass increases with r at the rate expected for a 



. Navarro, Frenk & White (1997) density profile. 

Or 



Subject headings: dark matter — galaxies: clusters: individual (Coma) — methods: 
\ statistical — survey 



•/"J ■ 1. INTRODUCTION 

For more than sixty years, the Coma cluster has been the most widely studied system of galaxies 



(see the review by Biviano 1998] ) . ^wicky (1933)| first estimated the mass of a system of galaxies by 



applying the virial theorem to Coma. This measurement provided the first piece of evidence that 
dark matter dominates the gravitational dynamics of systems on megaparsec scales. More accurate 
analyses use the Jeans equation to include a surface term in the virial theorem which accounts for 



the unsampled portion of the system (|The Sz White 1986|) . Major uncertainties in the mass arise 



from the assumption that light traces mass ( The fe White 19*8^ ; Merritt 1987] ) . Alternatively, X-ray 



observations provide information about the gravitational potential of the cluster and therefore, by 
assuming hydrostatic equilibrium, the total mass profile QHughes 1989| ; [Watt et al. 1992| ; |Hughes| 



1998). In spite of the assumptions involved, the mass determinations from optical and X-ray data 



in the core of the cluster (r <^ 0.5ft. 1 Mpc) agree within a factor of two. 

Although studies of the core of Coma are abundant, there are few quantitative studies of the 



outskirts (see e.g. the review by West 1998). These external regions are far from virial equilibrium 
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and therefore the data are more difficult to interpret. Furthermore, their large angular extent on 
the sky demands substantial amounts of telescope time for complete spectroscopic and photometric 
observations. Based on shallow (and often incomplete) redshift surveys, there have been analyses 
of this infall region to extract the value of the cosmological density parameter S7o ( [Bchectman 1982] : 
Capelato et al. 1982; Regos fc Geller 1989| ; van Haarlem et al. 1993). Alternatively, optical ( |Kent 



fc Gunn 1982| ; |Vedel fc Hartwick 1998|) and X-ray data ( Hughes 1989| ; Making 1994]) have been 
extrapolated to obtain a mass estimate at r ~ 3/i _1 Mpc. 

Galaxies falling into the potential well of a cluster populate a well-defined region in redshift 
space. Diaferio &: Geller (1997)| show that, in hierarchical clustering scenarios, the amplitude of the 
caustics which bound this region is proportional to the gravitational potential of the dark matter 
halo. In fact, because clusters accrete mass through the aggregation of smaller systems, random 
motions rather than spherical infall determine the amplitude of the caustics. Thus, measurement 
of the amplitude of the caustics yields the system mass, as long as galaxies trace the velocity 
field of the infall region reliably. This argument applies at radii as large as ~ lO/i -1 Mpc. N- 
body simulations of Cold Dark Matter models which include a semi-analytic treatment of galaxy 
formation and evolution suggest that there is no velocity bias on the relevant scales ( [Kauffmann 
et al. 1999| ; piaferio et al. 1999j ). piaferio (1999) uses these simulations to show that the mass 
estimation technique we apply here yields estimates to r = 5 — lOh^ 1 Mpc with a typical uncertainty 
of 50% for samples of a few hundred galaxies. 

As a first application of the mass estimation technique, we use a new redshift survey covering 
the Coma cluster to a 10° radius to measure the cluster mass profile. The largest of our sam- 
ples contains 1693 redshifts. We briefly describe the redshift survey (Sect. ||), outline the mass 
estimation technique (Sect. ||), and discuss the results (Sect. Q). 



2. THE DATA 

To explore the infall region surrounding the Coma cluster we measured 447 new redshifts for 
galaxies with m^ < 15.4 within 4.25° of the center of the Coma cluster and in a strip between 
26.5° < 5i95o < 32.5°. Our redshift survey is complete to niR = 15.4 within 4.25° of the cluster 
center. Within 10° of the cluster center the survey is complete to iur = 15.4 for 42% of the area 
and to mzw = 15.5 for the remaining area. 

We measured redshifts with the FAST spectrograph on the 1.5-meter telescope at Whipple 
Observatory. The FAST observations are typically 5-20 minute exposures with a 300 line mm -1 
grating, providing spectral coverage from 3600-7600A at 6A resolution. We use the IRAF task 
XCSAO in the RVSAO package (Kurtz and Mink 1998) to determine radial velocities based on 
absorption- and emission-line template cross-correlation. We compute the uncertainties in the 
velocities according to the procedures in Kurtz and Mink (1998). 

We identified galaxies brighter than mjj = 15.4 from Space Telescope Institute scans of the 
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POSS E plates (provided by B. McLean and M. Postman). We calibrated the instrumental magni- 
tudes by comparison with the Century Survey (Geller et al. 1997). For each plate in the strip (5 
plates) we fit the zero point offset and the non-linear scale error. For plate E1393, which contains 
the cluster center, we also calibrated against the CCD photometry of Kashikawa et al. (1998). For 
the Century Survey calibrators we obtained a zero point offset (at R = 15.0) of —0.13 mag and a 
slope for the scale error of —0.228 mag/mag; the Kashikawa calibrators yielded —0.16 and —0.223, 
respectively. We estimate that the magnitude limit is uncertain by ~ 0.1 m . The catalog to a radius 
of 10° contains all of the galaxies in the Zwicky (Zwicky & Herzog 1963) catalog with mz w < 15.5. 

We consider three different samples: (1) L4.25: 485 galaxies brighter than m^ = 15.4 within 
a 4.25° radius; (2) L10.0: 1077 galaxies with tcir < 15.4 or mz w < 15.5 within a 10° radius; and 
(3) C10.0: a compilation of 1693 galaxies within 10.0° regardless of their apparent magnitude. The 
C10.0 sample includes redshifts from Thorstensen et al. (1989), Wegner et al. (1990), Huchra et 
al. (1992), Thorstensen et al. (1995), van Haarlem et al. (1993), Colless & Dunn (1996), Geller et 
al. (1997), Falco et al. (1999). 



3. MASS PROFILE 

In redshift space, galaxies around clusters appear within regions delimited by caustics with 
a characteristic trumpet shape (Kaiser 1987): outside the caustics, the density of galaxies drops 



substantially. Galaxies outside the caustics are background or foreground galaxies (e.g. den Hartog 
k Katgert 1991 ). 



Half of the redshift space distance between the upper and the lower caustic at projected 
separation r from the cluster center defines the amplitude A(r). If we assume spherical symmetry, 
A(r) is a measure of the gravitational potential (j){r). We can thus estimate the cluster mass profile 
from the inner halo to the outer infall regions. In spite of the assumption of spherical symmetry, 
the technique yields mass profiles accurate to ~50% for a suite of iV-body simulations (Diaferio & 
Geller 1997; Diaferio 1999). 

We now apply the mass estimation technique described in Diaferio (1999) to Coma. To derive 
a cluster center and velocity, a hierarchical cluster analysis identifies the cluster members. An 
adaptive kernel method locates the peak of the density distribution of the cluster members on the 
plane of the sky; the cluster redshift is the median of the cluster member velocity distribution. 
In the cleanly magnitude limited sample L4.25, we identify 263 cluster members. The center we 
derive for Coma is a c (2000) = 12 h 59 m 25.0 s , 5 C (2000) = 27°56'45.7", and cz c = 7090 km s _1 . This 
position differs by a few arcminutes and by ~ 250 km s _1 from, for example, the center of the main 
condensation in Coma identified by Colless and Dunn (1996). The meaning of our center is not quite 
the same as earlier optical and X-ray determinations; ours is an estimate of the minimum of the 
global Coma system potential well derived from a more extended survey than available heretofore. 
We use the same center for all of the samples we analyze. 



- 4 - 



At the threshold in the binary tree used to identify the cluster members, there are a number 
of groups distinct from the main cluster. We use the procedure applied to the whole cluster to 
determine the center of each group. We then limit our analysis to individual galaxies, i.e. galaxies 
that do not belong to any group, and to galaxies in groups which lie within cz max — 4000 km s" 1 
of the cluster center. Larger values of cz max do not change our results; galaxies at these redshift 
distances have no effect on the location of the caustics. Smaller values might exclude galaxies in 
the cluster core. We locate all galaxies in the redshift diagram (r, v), where 

r = — -sin6>; v = cz — cz c cos 8. (1) 
Ho 

Here, cz c , cz, and 6 are the redshift of the cluster center, the galaxy redshift and the angular 
separation between the cluster center and the galaxy, respectively. 

We now use an adaptive kernel method to compute the two-dimensional density distribution 
f q (r, v) of the galaxies within the redshift diagram. To compute f q (r, v) with a spherical smoothing 
window, we rescale r and v so that the ratio h v /h r of the smoothing window sizes along v and r 
respectively is q = 25. This choice yields equal weight to the typical uncertainties on v and r, e.g. 
50 km s _1 and 0.02/i _1 Mpc for nearby clusters, respectively. Different values of q in the range 
[10, 50] have negligible effects on the results. We then identify the threshold k which minimizes the 
quantity 

S( K ,R) = \(vl c ) K , R -4(v 2 ) R \\ (2) 

where {v 2 )r is the velocity dispersion of the cluster members and {v 2 sc ) k ,r is the escape velocity 
computed with the estimated A(r); R = l.Wh^ 1 Mpc is the mean distance of the cluster members 
from the cluster center. We compute the caustic amplitude A(r) = min{t> u (r), Vd(r)} where v u (r) 
and Vd(r) are the solutions of the equation f q (r, v) = k. 

The mass M(< r) within the radius r from the cluster center is 

GM(< r) = l r A 2 (x)dx. (3) 
2 Jo 

To quantify the error in the amplitude and the mass profile, we assume that the relative error 
5A(r)/ A(r) = k/ max{/ g (r, v )} and <5Mj = J2j=i,i \^ m j^A(rj)/A(rj)\, where mj is the mass of the 
shell [rj—i, rj] 

Gmj = - f 1 A 2 {x)dx. (4) 

A^-body simulations show that this expression approximates the 1-a spread of profiles obtained by 
projecting the clusters along different lines of sight. 

The upper panels of Figure |l] show the redshift diagrams of Coma for our three samples. Each 
case shows a clearly defined high density region around the median redshift of Coma. The bold 
lines show where our method locates the caustics. Variations of A(r) from sample to sample are 
within the 3-cr uncertainty (middle panels). Thus, the mass profile (eq. H) is also consistent from 
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sample to sample. The shaded areas in the middle and bottom panels show the 2-a uncertainty of 
the profiles. Because our three samples contain galaxies of different luminosity, the consistency of 
the profiles indicates that (1) galaxies are tracers and that (2) velocity bias is negligible as we have 
assumed. 

We limit our profile to r = lO/i -1 Mpc. At larger distances, the meaning of a "cluster" 
is unclear. Our largest sample indicates M(< 5.5/i _1 Mpc) = (1.44 ± 0.29) x 10 15 /i _1 M© and 
M(< 10/i _1 Mpc) = (1.65 ± 0.41) x 1O 15 /i _1 M (1-a error). Error bars in Figure 1 show the range 
of X-ray mass estimates available at r = 0.5/i _1 Mpc and r = 2.5h~ 1 Mpc ( Hughes 1989| ). At these 
small radii, our estimate agrees very well with the X-ray estimates. 

We also plot the cumulative mass profile for a softened isothermal sphere (short-dashed line) 
and for a halo with a Navarro, Frenk & White (1997; NFW) profile (long-dashed line). We fit these 
respective mass profiles M(< r) = 4npor^[(r/r c ) — arctan(r/r c )] and M(< r) = 47r<5 c/ o c rf {log[l + 
(r/r s )] — (r /r s )/[l + (r /r s )]} in the range r = [0, l]/i _1 Mpc; po, r c , 5 C and r s are fitting parameters, 
and p c = 3Hq/8ttG is the critical density of the Universe. For the NFW profile we find r s = 
0.182 ± 0.030, 0.167 ± 0.029, and 0.192 ± 0.035/j" 1 Mpc for the L4.25, L10.0, and C10.0 sample, 
respectively. Note that larger fit ranges, up to r = [0,5]/i _1 Mpc, do not change r s . The softened 
isothermal sphere yields a poor fit for any fit range larger than r = [0, 2]/i _1 Mpc. The fit parameters 
determine the expected behavior of the profiles at larger radii. Our measurement shows clearly that 
the NFW extrapolation agrees remarkably well with the observed mass increase. The mass profile 
for the isothermal sphere increases too steeply for consistency with the data. 

The parameter 5 C in the NFW profile is related to the concentration parameter c = r2oo/r s , 
where r2oo is the radius of the sphere with average mass density 200 times the critical density. We 
can therefore estimate the virial radius r2oo of Coma: we find r2oo — 1-5/t Mpc for each of our 
three samples. 



4. DISCUSSION 

We use a redshift survey within 10° of the center of the Coma cluster to measure the mass 
profile of the cluster to r = lO/i" 1 Mpc: M(< lO/T^Mpc) = (1.65 ±0.41) x 10 15 /i _1 M o (1-a error). 
This mass profile extends to distances 50 times larger than the cluster core radius r s ~ 0.2/i -1 
Mpc and to substantially larger radius than previous estimates. The mass profile appears robust: 
it does not change appreciably when we increase the size of the dataset. 

The mass profile increases at a rate consistent with the NFW density profile and substantially 
smaller than the density profile for an isothermal sphere. Because the NFW mass profile increases 
only logarithmically, our results do not exclude convergence of the mass of Coma at r ;> 3 — 4/i _1 
Mpc. However, the agreement between our measurement and the increasing NFW profile does 
suggest that, other than a density contrast criterion, there is no obvious definition of the extent of 
the halo of Coma. We can use the NFW density profile fit to estimate the virial radius of Coma 
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1.5/1" 1 Mpc. 



Geller et aL (1999) will include the data, an analysis of the substructure within the infall region 
of Coma, and a discussion of the morphologies and spectroscopic properties of galaxies in the infall 
region. In forthcoming papers, we will apply our mass estimation technique to similar data obtained 
for the infall regions of the Abell clusters A539, A576 and A1367 and for the A2199/A2197 double 
cluster. 

The mass estimation method we apply here depends on kinematic data only. If velocity 
biases are weak, as suggested by high resolution A-body simulations including a phenomenological 



treatment of galaxy formation ( Kauffmann et al. 199E ; Diaferio et al. 1999 ), the mass estimate is 
independent of the luminosity distribution: with accurate photometric data we can then measure 
the behavior of the mass-to-light ratio out to ~ 5 — lO/i^ 1 Mpc scale. In the future we hope to 
report such a measurement for the Coma cluster. 

From photometric data, we can also quantify the variation of the galaxy population and its 
kinematic properties with the distance from the cluster center. Within the virial radius of clusters 



the fraction of blue galaxies increases with radius (Dressier 1980; Whitmore & Gilmore 1992 



Carlberg et al. 1997| ; |de Theije &: Katgert 1999 ) and their velocity dispersion is generally larger 



than the velocity dispersion of red galaxies ( Mohr et al. 1996| ; |Carlberg et al. 1997 ), indicating 
that blue galaxies are actually falling into the cluster and do not satisfy virial equilibrium (e.g. 
Biviano et al. 19971 ). A-body simulations of Cold Dark Matter models naturally reproduce these 



properties ( piaferio et al. 1999 ). These simulations also show that kinematic differences between 
blue and red galaxies become smaller at radii larger than the virial radius, but at different rates in 
universes with different £Iq. Therefore, by combining photometric and spectroscopic data covering 
the infall regions of galaxy clusters over a large redshift range we can potentially discriminate 
among cosmological models. 
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Fig. 1. — Top panels: Galaxy distribution in the redshift diagram for our three samples. The bold 
lines indicate the location of the caustics. There are 332, 480, and 691 galaxies within the caustics 
in the samples L4.25, L10.0, and C10.0, respectively. Half the distance between the caustics defines 
the amplitude A(r) shown in the middle panels. Bottom panels: The bold lines are the mass 
profiles derived from the data with eq. (|3|). The two error bars show the range of the X-ray mass 
estimates listed in Hughes (1989). Short-dashed and long-dashed lines are the cumulative mass 
profile for a softened isothermal sphere and an NFW density profile with parameters obtained by 
fitting the mass profile in the range [0, l]/i _1 Mpc. Shaded areas in the middle and bottom panels 
indicate the 2-a uncertainty. 



